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ABSTRACT : In this paper presents, a tuned parameters PID controller is designed for the original 
higher order system through its lower order reduced model. The PID parameters are obtained by using 
classical control technique and these parameters are tuned through Ziegler’s Nichols method using 
MATLAB software. The designing of PID controller can improve settling time and steady state response. 
Designing of a controller or compensator for the original higher order system is complex and costly. So to 
avoid the complexity, first controller is designed for the reduced order model of the original system and 
then this designed controller is connected to the original system. Defiantly, this Designed PID controller 
can improve the time response characteristics of the original higher order system. So in this paper for 
determining the reduced order model, a mixed method is used. The combination of PSO and Polynomial 
are the methods used to reduce the high order system to low order system. PSO technique is used for 
determining the denominator coefficients and Polynomial method is used for obtaining the numerator 
coefficients. The original and reduced order systems without PID controller and with PID controller are 
tested with numerical examples. 

Key Words: PID controller, Ziegler’s Nichols method, PSO optimization, polynomial technique, order 
reduction, transfer function, settling time, steady state error. 



I. Introduction 

The order reduction techniques play a very important role in the control systems. The design of controllers and 
compensator for original system is costly and tedious. . In order to perform simulation analysis or control design 
on those higher order models one will face many difficulties. Designing of compensator and controller to the 
reduced order system is easy and simple. So to avoid the above problem order reduction is necessary. 

There are many order reduction methods proposed in both time domain and frequency domain for 
reducing order from the higher order system to lower order system. One of the order reduction method is pade 
approximation, Shamash [1] proposed this method. The disadvantage of the pade method is that once the 
reduced order model is unstable, if the poles of the original higher order systems are more, then reduced order 
model is checked the stability for every time. Another method is Routh approximation, Hutton and Friendland 
[2] proposed this method. The drawback of the rough approximation method is no guarantee of the stability of 
the reduced order system. 

To reduce the order from higher order to lower order system Evolutionary techniques [6] are used. So 
many Evolutionary techniques are there, among all the Evolutionary techniques Particle swarm optimization 
technique is most useful technique. Particle swarm optimization technique is similar to the Genetic algorithm 
technique. The difference between the PSO and Genetic algorithm is the PSO is an easy and simple algorithm. 
In this paper a tuned parameter PID controller is designed for the original higher order system through its lower 
order system. By designing the PID controller, we have to improve the settling time and steady state error of the 
original and reduced order systems. The reduce order model is obtained the mixed methods. The mixed methods 
are a PSO and Polynomial method. 



II. PROCEDURE FOR PSO AND POLYNOMIAL METHODS 

II. 1. Procedure for PSO: 

The original higher order system transfer function of order ‘n’ is 

C (S) W< ~ S) ao+ais+a 2 s 2 +- a m _ 1 5"~ 1 

n D(s) />0+/)lS+/)2S 2 +-" t) n S n 

The reduced order system transfer function of order ‘k’ is 

R (o') = N AO d 0 +d 1 s+d 2 s 1 + • • i t k 1 

k ^ BfcO) e 0 +e 1 s+e 2 s 2 +-- e k s k 

The PSO method is the one very important optimization method for solving the optimization problems. It is the 
swarm intelligence technique and computer oriented algorithm. The particles are the very important for the PSO 
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algorithm, these particles are randomly generated from initial population. These particles are moving in the 
search space according to updating its own position and velocity. Every particle in the algorithm is affected by 
its personal best (pbest) that is also called as local best and global best (g best) [13] for obtaining the best 
position in the search space. In this PSO algorithm for every iteration the values of personal best and global best 
are calculated and by using these values position and velocity of the particles are updated. For each iteration the 
best position values are called personal best (pbest) and among all the iteration the best fitness value in the 
swarm or population is called global best (gbest).the velocity and position values are updated by using below 
formulae. 



v id +1 = v id + c i * r i * (PUL - x-d) + c 2 *r 2 * ( p k gd - x k d ) ( 3 ) 

x?d +1 = 4 + 4 (4) 

Where Ci and c 2 are the cognition and social parameters and i'i and r 2 are the random parameters. Ci and c 2 
values generally taken as ‘2’ these parameters are either variable or constant. The ranges of random parameters 
are in between 0 and 1. And p id and p gd are personal best and global best values. For balancing the personal best 
or local best and global best values inertia weight ‘w’ is introduced. 

w = Wf + {wf — Wj)(maxtt — it) /it (5) 

Then the above equation 3 can be written as 

Vut 1 = w * v k d +c 1 *r 1 * (p'/ d - x? d ) + c 2 *r 2 * ( p k gd - x* d ) (6) 

11.2. Procedure for Polynomial method: 

For finding the values of the numerator coefficients polynomial method is used. The numerator coefficients are 
obtained by equating the original higher order system and reduce order system. 

ao+ais+a2S 2 +-" a n _is n_1 _ do+dis+d 2 s 2 +'" d^-js^ -1 

fao+Z>ls+/)2 s2 +’" ti n s n eo+eis+e2S 2 +--- e k sk 



Then equate the power’s of ‘s’ on both sides 

a 0 e 0 = ^0^o (8) 

a 0 e 1 + ct-^e q= hod 2 + h 2 dg (9) 

ci 0^2 + ct\6\ T a 2 e 0 = bQd 2 T bid'ib 2 dQ (10) 

a n-l e k~b n dk~l (11) 



By solving the above equations we get the values of d 0 ,dj and d 2 ,....d q 

III. PID Controller 

The combinations of suitable parameters that are proportional, integral and derivative controller can improve the 
all aspects of system performances. The prortional controller settle the gain, but produces a steady state error. 
The integral controller decreases or avoid the steady state error The derivative controller reduces the rate of 
error. By using this PID controller we can improve the settling time and steady state errors. 




Fig-1: PID controller Block diagram 

In order to obtain the output of the above controller, the P, I and D controller terms are added. Then the final 
expression of PID controller is, 

u(t) = K p e(t) + K t J 0 f e(T)dT + K d |e(t) (12) 

In this paper, PID controller is designed and tuned by using Ziegler’s Nichols method. The PID controller 
parameters are tuned by classical control technique through Ziegler’s Nichols method in MATLAB software. 
The tuned parameters of PID controller using Ziegler’s Nichols method are shown in below table. 
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Table 1: Values P, PI, PIP using Ziegler’s Nichols method 



Type of Controller 


Kp 


Ki 


K d 


P 


0.5K„ 


00 


0 


PI 


0.45K U 


1.2 K p 
T„ 


0 


PID 


0.60K,, 


2K p 

L 


KpT u 

8 



IV. Numerical Examples 

IV. 1. Example: 

IV. 1.1. Determination of reduced order system: 

Consider the original higher order transfer function of an order ‘5’ [6] 

. 10s 4 +82s 3 +264s 2 +396s+156 

G(s)= 



(13) 



2s 5 +21s 4 +84s 3 +173s 2 +148s+40 

Numerator coefficients are obtained by PSO algorithm. For implementing this several parameters are to be 
considered. The cognition Ci and social parameters c 2 are ‘2’ and the range of random parameters are in 
between 0 and 1 . 

Population size or Swarm size = 30 
Coefficients to be finding = 2 
Maximum number of iterations = 300 

The obtained denominator polynomial for second order reduced model by using PSO is 

D 2 (s) = 2.36919 + 3.918816s + s 2 (14) 

For finding the values of the numerator coefficients polynomial method is used. The numerator coefficients are 
obtained by equating the original higher order system and reduce order system 

10s 4 +82s 3 +264s 2 +396s+156 rfo+rf^s ^ 

2s 5 +21s 4 +84s 3 +173s 2 + 148s+40 _ 2.36919+3.91886s+s 2 ' ’ 

On equating the power’s of ‘s’ on both sides we get numerator coefficients and multiplying numerator with k to 
match the steady state response . 

Therefore, the numerator polynomial by polynomial method is 

N 2 (s) = 9. 255+5. 717s (16) 

The obtained reduced order model by using PSO and Polynomial method is shown below. 

. . 9.255+5.717s 

R+s) = 7 (17) 

“ v ' 2.369 19 +3.9 1886 s+s 2 v ' 



TV. 1.2. PID Controller design using reduced order model: 

For designing PID controller first reduced original higher order system into lower order system. The obtained 
reduced order model by using PSO and Polynomial method is shown below. 

9. 255+5. 717 s 

R+S = 7 (18) 

J 2.36919 +3.91886s+s 2 v ' 

For the above reduced order model the PID controller is designed. The tuning parameters of PID controller are 
obtained by using Ziegler’s Nichols method: 

K p = 2.24096 K, = 14.006 K d = 0 
The transfer function of PID controller is 

G d O) = (19) 

The reduced order system with PID controller is 

. , 12.81 s 2 +100.8s+129.6 

R-> r (s) = -5 5 (20) 

' s 3 +16.73s 2 +103.2s+129.6 v 

The original high order system with PID controller is 

22.41 s 5 +323.8s 4 +1740 s 3 +4585 s 2 +5895 s+2 185 

2s 6 +43.41s 5 +407.8s 4 +1913s 3 +4733s 2 +5936s+2185 ( ’ 



Table 2: comparison of time respons characteristics with and without PID controller 





Settling time (Sec) 


Rise time (Sec) 


Steady state value 


Original higher order system 
without PID 


5.45 


2.76 


3.9 


Original order system with PID 
controller 


0.68 


0.124 


1 


Reduced order system with PID 


4.77 


2.54 


3.91 


Reduced order system with PID 
controller 


0.616 


0.112 


1 
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System: Reduced w ithout PID 
4 Settling Time (sec): 4.77 



Step Response 



System: Reduced w ith HD 
Settling Time (sec): 0.616 






System: Reduced w ithout HD 
Final Value: 3.91 



System: Reduced w ith HD 
Final Value: 1 



- Reduced w ithout PID 
Reduced w ith PID 



Time (sec) 



Fig-2: Step responses of reduced order model with and without PID controller 



Step Response 



System: Original w ithout HD 
Rnal Value: 3.9 



System: Original w ithout HD 
Settling Time (sec): 5.45 



System: Original w ith PID 
Settling Time (sec): 0.68 



System: Original w ith HD 
Rnal Value: 1 



- Original w ithout HD 

- Original w ith PID 



3 4 5 

Time (sec) 



Fig-3: Step responses of original order model with and without PID controller. 



IV.2. Example: 

IV. 2.1. Determination of reduced order system: 

Consider the original higher order transfer function of an order ‘6’ [4] 

. 3s 5 +108 s 4 +122 Is 3 +5792 s 2 +l 1860 s +8400 

G(s) = 



( 22 ) 



s 6 +41 s 5 +571 s 4 +349 Is 3 +10060 s 2 +13100 s+6000 

Numerator coefficients are obtained by PSO algorithm. For implementing this several parameters are to be 
considered. The cognition Ci and social parameters C 2 are ‘2’ and the range of random parameters are in 
between 0 and 1 . 

Population size or Swarm size = 30 
Coefficients to be finding = 2 
Maximum number of iterations = 300 

The obtained denominator polynomial for second order reduced model by using PSO is 

D 2 (s) = 8.714992 + 8.630098s + s 2 (23) 

For finding the values of the numerator coefficients polynomial method is used. The numerator coefficients are 
obtained by equating the original higher order system and reduce order system 

3s 5 +108 s 4 +122 Is 3 +5792 s 2 +11860 s+8400 _ d 0 +dis 

s 6 +41 s 5 +571 s 4 +349 Is 3 +10060 s 2 +13100 s+6000 “ 8.714992 +8.630098 +s 2 

On equating the power’s of ‘s’ on both sides we get numerator coefficients and multiplying numerator with k to 
match the steady state response 

Therefore, the numerator polynomial by polynomial method is 
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N 2 (s) = 12.2039 + 3.355s (25) 

The obtained reduced order model by using PSO and Polynomial method is shown below. 

n f >. 12. 2039+3. 355s 

Ry is) — (26) 

z w 8.714992 +8.630098 s+s 2 v ' 



IV.2.2. PID Controller design using reduced order model: 

For designing PID controller first reduced original higher order system into lower order system. The obtained 
reduced order model by using PSO and Polynomial method is shown below. 

3.355s+12.2039 x 

R^Cs) — — 5 (27) 

2V ’ s 2 +8.630098s+8.714992 v ’ 

For the above reduced order model the PID controller is designed. The tuning parameters of PID controller are 
obtained by using Ziegler’s Nichols method: 

K p = 6.46656 K, = 53.888 K d = 0 
The transfer function of PID controller is 

53.888+6.46656s /oo\ 

— : (2o) 



Gd( s ) — 

The reduced order system with PID controller is 

. . 21.7s 2 +259.75s+657.6 

R_2c(s) — 9 

2 s 3 +30.335 s 2 +268. 4s +657.6 

The original high order system with PID controller is 

19.4s 6 +860.1s 5 +13720 s 4 +103300 s 3 +3888000 s 2 +693400 s+ 452700 
_ s 7 +60. 4s 6 +1431 s 5 +17210 s 4 +113300 +401900 s 2 +699400 s+452700 



(29) 



(30) 



Table 2: comparison of time response characteristics with and without PID controller 





Settling time (Sec) 


Rise time (Sec) 


Steady state value 


Original higher order system 
without PID 


3.47 


1.84 


1.4 


Original order system with PID 
controller 


0.529 


0.0839 


1 


Reduced order system with PID 


3.16 


1.74 


1.4 


Reduced order system with PID 
controller 


0.517 


0.0761 


1 



1.4 



1.2 



1 

<d 0.8 

"O 
“O- 

E 

< 0.6 
0.4 
0.2 
0 

0 1 2 3 4 5 6 

Time (sec) 

Fig-4: Step responses of reduced order model with and without PID controller. 



Step Response 



System: Reduced w ithout PID 
Settling Time (sec): 3.16 



System: Reduced w ithout PID 
Final Value: 1 .4 



System: Reduced w ith PID 
Settling Time (sec): 0.517 



System: Reduced w ith PID 
Final Value: 1 



Reduced w ithout PID 
Reduced w ith PID 
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Time (sec) 



Fig-5: Step responses of original order model with and without PID controller 

V. Conclusion 

In this paper, a tuned parameters of PID controller is designed for original higher order system through 
its lower order reduced model. The PID parameters are obtained by using classical control technique and these 
parameters are tuned through Ziegler’s Nichols method using MATLAB software. The designed PID controller 
is connected to original and reduced order systems. The design of PID controller improves the settling time and 
steady state error than the proposed reduced order method. The combination of PSO and Polynomial are the 
methods used to reduce the high order system to low order system. PSO technique is used for determining the 
denominator coefficients and Polynomial method is used for obtaining the numerator coefficients. The original 
and reduced order systems without PID controller and with PID controller are tested with numerical examples. 
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